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Abstract 

The equations for geodesic flow on the ehipsoid are weh known, and were first 
solved by Jacobi in 1838 by separating the variables of the Hamilton-Jacobi equa- 
tion. In 1979 Moser investigated the case of the general ellipsoid with distinct 
semi-axes and described a set of integrals which weren't know classically. After 
reviewing the properties of geodesic fiow on the three dimensional ellipsoid with 
distinct semi-axes, we investigate the three dimensional ellipsoid with the two mid- 
dle semi-axes being equal, corresponding to a Hamiltonian invariant under rotations. 
The system is Liouville-integrable and thus the invariant manifolds corresponding 
to regular points of the energy momentum map are 3-dimensional tori. An analysis 
of the critical points of the energy momentum maps gives the bifurcation diagram. 
We find the fibres of the critical values of the energy momentum map, and carry out 
an analysis of the action variables. We show that the obstruction to the existence 
of single valued globally smooth action variables is monodromy. 



1 Introduction 

The geodesic flow on the ellipsoid is the classical example of a non-trivial separable and 
thus Liouville integrable Hamiltonian system. It is the prime example in Jacobi's famous 
"Vorlesungen iiber Dynamik" |H] and may be considered as his motivation to develop 
Hamilton-Jacobi theory and the solution of the Abel- Jacobi inversion problem. Its modern 
treatment was pioneered by the Zurich school, namely by Moser jT2| and Knorrer llOj. 
generalising to the n-ellipsoid and providing smooth integrals and the general solution in 
terms of ^-functions for the generic case of an n-ellipsoid with pair-wise distinct semi- 
axes. Separation leads to a curve of genus n and the n actions are given by integrating a 
differential of second kind over a basis of real cycles. The generic motion on an ra-torus 
corresponds to a non-degenerate curve. Special motions correspond to degenerate curves. 

*einail: c.ni.davison@lboro.ac.uk, h.r.dullin@lboro.ac.uk, a.bolsinov@lboro.ac.uk 
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For two degrees of freedom the topology of the Liouville fohation was analysed in pp. 
An excellent general approach to the topology of Stackel systems, including geodesic flow 
on the ellipsoid with distinct semi-axes, was carried out by Zung This is the first 
paper where the singularities of the Liouville foliation were studied for geodesic flow on 
the ellipsoid. We extend his results to the degenerate case. 

Surprisingly most of these results are not stable when the ellipsoid becomes degenerate, 
i.e. when some semi-axes coincide. The smooth integrals develop poles in this limit, the 
hyperelliptic curve changes, complete Abelian actions integrals change from 2nd kind to 
3rd kind, and the topology of the critical values in the image of the energy-momentum 
map changes. Here we study these changes for the 3-ellipsoid. The most interesting result 
appears when the middle-axes coincide; the set of regular values of the energy-momentum 
map becomes non-simply connected. Duistermaat j3] realised that in this case global 
action variables might not exist. The torus-bundle over the regular values is non-trivial 
and has monodromy. This shows that by making the system simpler (i.e. more symmetric) 
it can become more complicated (i.e. have a non-trivial torus bundle). 



2 The geodesic flow on generic S-eUipsoids 

A 3-ellipsoid embedded in M"^ with coordinates x = (xq, Xi, 0:2, Xs) has the equation 
{A~^x,x) = 1 with the standard Euclidean scalar product (,) and a positive definite 
matrix A. This quadratic form can always be diagonalised by an orthogonal transfor- 
mation and the eigenvalues of A are denoted by < ao < «i < «2 < as- Thus we 
may assume that A is diagonal and the standard form of the ellipsoid with semi-axis ^/oi 
embedded in is 

2222 

/-y»ij rf'-^ rr>'-' rf'^ 

ao a\ a2 0:3 

For the generic non-degenerate ellipsoid the semi-axes are distinct. The Lagrangian of 
a free particle with mass 1 is L = j{xq + x1 + x^ + x|). The equations of motion with 
Lagrange multiplier A enforcing the constraint Ci = are 

X = ^AA-'x, A = . j^''^' ^] . . (1) 

{A'^x,A-^x) ^ ' 

A Hamiltonian description can be obtained by introducing momenta yi = Xi and enforcing 
the constraint by replacing the standard symplectic structure dx Ad|/ by a Dirac bracket. 
The Dirac bracket has as Casimirs the constraint for being on the ellipsoid Ci = and 
the constraint for its tangent space 

_ xoyo xiyi X2y2 , xsys _ 

O2 — 1 1 1 — u . 

ao cti <^2 tts 

With the notation 

^ Xn X? Xa Xo 1 dCi dC2 

D = — + — + — + — = -y — - — - 

c^o al al al 2 ^ dxi dyi 
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the Dirac bracket with Casimirs Ci and C2 is given by 

{a;.,x,} = 0, {a:„,,} = 5.,--i^, {^^^ = "^^^f^ • (2) 

The Hamihonian is H = ^{Vq + yf + 1/2 + vl) and the equations of motion are 

Xi = {xi,H}, yi = {yi,H}, 2 = 0,1,2,3. (3) 

These equations are equivalent to (^. The Hamiltonian vector field generated by H is 
denoted by Xh- 

The system is Liouville integrable with smooth global integrals (in the generic case of 
distinct semi-axes) first found by Uhlenbeck and Moser jT2] 

F.=y^+ Yl ' > . = 0,...,3. (4) 

On the symplectic leaf of the Dirac bracket given by Ci = C2 = they are related by 
Y^Fi/ai = and they have pair- wise vanishing brackets jl2] . The integrals Fi are related 
to the Hamiltonian by = i(Fo + Fi + F2 + F3). 

Ellipsoidal coordinates are local coordinates on the ellipsoid that separate the Hamil- 
tonian. They are defined as the roots A of K{x, x; X) = 1 where 



K{x,y;X) = J2 



•^iyi 



at - X 

The equations K{x, x; Aj) = 1 are linear in xf and can be easily solved to give 
2 B{ai) 



X, 



B{z) = ll{X,-z), A{z) = ll{a,-z). (5) 



A'iat) 

Because of the poles in K for fixed x the 4 roots satisfy 

Aq ^ c^o ^ -^1 ^ ^ -^2 ^ ^2 ^ A3 < 0^3 . 

Fixing Ao = gives a coordinate system (Ai, A2, A3) on the ellipsoid since K{x, x; 0) — 1 = 
Ci- The coordinate transformation to the new variables Aj and their conjugate momenta 
Pi gives 

V2^^n;(A,-A.)- 

The primed product excludes the vanishing term with j = i. The geodesic fiow on the 
ellipsoid Aq = is described by the invariant subset given by po = Aq = 0. The variables 
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can be separated by using the van der Monde matrix (A^ )jj as a Stackel matrix [Zj. 
With separation constants Sj where S3 = 2h and Sq = the separated equations are 



P-=-^^, Q{z)=2hz-' + S2z' + siz + s^. (6) 

The system separates on the hypereUiptic curve w"^ = —Q{z)A{z) of degree 7, hence genus 
3. The relation between Fj and the separation constants Sj is determined by the residues 
of the identity 



' Fi Q{z) 



i=0 



z — ai A{z) 



(7) 



In particular Sa = - Xlcyd /o(«i + «2 + 03), Si = Xlcyd /o(«ia2 + OiOs + 0203), where fi 
denotes a value of Fi. 

The ellipsoidal coordinates \j only determine the squares of the Xj and thus have 
singularities when = 0. Smoother coordinates 0j on a covering torus designed so that 
their cotangent lift cancels the singularities given by A{z) in (jHl) are defined by 

^ = —j= , V . / X N ' "i-i<^i<ai, i = l,2,3. (8) 
dAi 2^{-lYA{\) 

This defines elliptic functions \i{4>i) with modulus k"^ = (^3 — «2)(ai — ao)/(('^3 ~ 
cn)(a2 — cto)) given by the cross ratios of the semi-axes squared for i = 1,3 and with 
complementary modulus ^/l — k"^ for i = 2. The momenta conjugate to (pi are denoted 
by pi. In this coordinate system the squares of the new momenta are smooth functions 
p'^ = {—iy~^^Q{Xi{(j)i)). The Hamiltonian H and the constants of motions from separation 
5*1, 5*2 in these coordinates are 

"2 ^2 ^2 

^ ^ Pi ^ P2 , Pi 



2Ai(A2 - Ai)(A3 - Ai) 2A2(A2 - Ai)(A3 - A2) 2A3(A3 - Ai)(A3 - A2) ' 

_ X2hpl , A1A2P2 , A1A2P3 

Oi — - — T- — -7- — r + - — — — — r- — - + 



^2 



Ai(A2-Ai)(A3-Ai) A2(A2-Ai)(A3-A2) A3(A3 - Ai)(A3 - A2) ' 
(A2 + X3)pi (Ai + X3)pI (Ai + X2)pI 



Ai(A2 - Ai)(A3 - Ai) A2(A2 - Ai)(A3 - A2) A3(A3 - Ai)(A3 - A2)' 



In these formulas each Aj represents the elliptic function Aj(0i). 

The {4>,p) coordinate system still is not a global coordinate system on the cotangent 
bundle of the 3-ellipsoid (such global coordinates do not exist). It has singularities at 
the "umbilical curve" ^ determined by Ai = A2 = ai contained in the Xi = plane and 
A2 = A3 = 02 contained in the X2 = plane. Explicit formulas are obtained by inserting 



^The term "umbilical curve" is used in analogy to the umbillic points on the two dimensional ellipsoid 
we don't know whether it has a differential-geometric characterisation. 
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Figure 1: Bifurcation Diagram of the generic Ellipsoid with = 1/3, 1,3,4 

these conditions into where A3 or Ai become the curve parameter, respectively. Hence 
the umbilic curves are coordinate lines on the respective sub-ellipsoid. The umbilic curves 
are thus 2 topological circles in the X2-a;3-plane with xq > or xq < and 2 topological 
circles in the Xo-Xi-plane with X3 > or X3 < 0. At these points in configuration space 
H and Si are singular for arbitrary momenta and a different coordinate system needs to 
be used. 

There are four 2-dimensional totally geodesic subflows obtained from setting Xi = yi = 
0. Similarly there are six 1-dimensional totally geodesic subflows obtained from setting 
Xi = Xj = yi = yj = for each pair of indices with i < j. A 1-dimensional subfiow is 
of course simply two periodic orbits. These six times two periodic orbits and the four 
subflows on 2-ellipsoids give the backbone of the bifurcation diagram, i.e. the critical 
values in the image of the energy-momentum map. These critical values are shown in 
figure m 

Each separated equation — (— l)*+^(5(Ai(0i)) = defines a curve in the phase portrait 
in the {(pi, pi) plane. The critical points occur at ]3j = and (pi determined by Q'{\i)\[ = 0. 
Solutions are of two types: Either double roots of Q or critical points of the elliptic function 
Aj. According to (jHl) critical points of Aj(0j) occur exactly for A, = or Aj = ctj. By © 
this implies that Xi-i = or Xj = 0, respectively, and similarly for yi from the cotangent 
lift of (jHl). These critical points therefore correspond to the geodesic subflows mentioned 
above. The image of the critical points Xj = yj = is given by a segment of the line in 
the S1-S2 plane given by Q{aj) = 0, see figure [T] The other type of critical points occur 
for those values of Si,S2 for which there is a double root in Q{z) = 2hz{z — with 
«i < < ^2 so that A2 = (i is fixed for this motion. For d not at its boundary values 
these critical points are not contained in any geodesic subfiow. 

Finally we have to establish whether any point with x-coordinates in the umbilic curves 
is critical. The umbilic curves are contained in xi = and X2 = 0, respectively. The 
gradient of Fi vanishes on the subfiow Xi = yi = 0. This shows that points on the umbilic 
curve contained in Xj = which have vanishing momentum yi = are indeed critical. 
But what about other momenta? If the momentum yi is non-zero the corresponding 
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Figure 2: Bifurcation Diagram of almost degenerate Ellipsoids with = (1/4,1/4 + 
£, 1, 2), (1/4, 1/2, 1/2 + £, 2), (1/4, 1/2, 1, 1 + e), (1/4, 1/4 + £, 1, 1 + e), where e = 0.03. 

geodesic will leave the plane Xj = 0. As soon as it is outside the sub-ellipsoid ellipsoidal 
coordinates are regular, and thus the geodesic is non-critical since the only critical points 
outside sub-ellipsoids are tori with fixed A2 with ai < A2 < 02; but these tori have no 
point in common with the umbilic curve unless A2 = «i or a2- 

We have thus proved the well known result that the bifurcation diagram is obtained 
from collisions of roots of the hyperelliptic curve w"^ = —Q{z)A{z). Subflows correspond 
to Q{z) having a root that coincides with a root of A{z). Hence the four sub-ellipsoids 
are given by the lines Q{aj) = in the image of the energy momentum map {h,S2,Si). 
For geodesic flows the energy can be fixed to 1/2 without loss of generality, and thus the 
four lines in (si, S2) space are the straight lines a| -|- S2aj + si = 0, see figured These four 
lines intersect in six points (si, S2) = {ajCtk, —C(j — ctk), corresponding to periodic motion 
in the jfc-plane. The other curved line of the bifurcation diagram is given by double roots 
in Q{z) = z{z — dy such that (si,S2) = {(P,—2d) where ai < d = X2 < a2 attaching 
tangentially to the straight lines of intermediate slopes at the codimension two points 
{af, —2ai), i = 1,2. The four disjoint regions of regular values have 2 or 4 tori in their 
preimage. 

Finally it is important to establish the (no n-) degeneracy and type |2] of the singular 
points in phase space. The gradient of Fi vanishes in the plane Xi = yi = since each term 
is quadratic and contains the vanishing terms. For non- degeneracy the spectrum of the 
Jacobian matrix of the flow of Fi, which is denoted by DXp. needs to be computed. When 
restricted to the plane Xt = i/i = the only nonzero entries are in the Xi-yi sub-block. 
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This sub-block reads 



fi^h'l '^"^tr^'^l where Ux,y) = T^^. (9) 

Notice that Ki{y,y) never vanishes when Fq and F3 [i = 0,3) are considered. Otherwise 
not all terms of the same sign and Ki{y,y) can vanish. Since this matrix is traceless the 
square of the eigenvalues is given by the negative determinant. This condition needs to be 
evaluated on a point of the singular fibre in question. Consider, say, Fi. Then the point 
2^0 = 2:3 = is on all critical sets. Then 1/2 = since this point is at a maximum of X2 on 
the sub-ellipsoid xq = 0; moreover X2 = ±^/a2- At this point the diagonal terms vanish, 
and the eigenvalues vanish when Ki{y,y) vanishes. For Fi = this occurs for the special 
momentum when y^ = 2h{aQ — Oi) / (ao — 03), and this is the point of tangency of the line 
Fi = with the curve of double roots in Q in the bifurcation diagram. All other points are 
non-degenerate. Similarly for F2. This shows that there are two degenerate singularities 
in the geodesic flow on the non-degenerate n-dimensional ellipsoid when n = 3. This is 
why the claim made in Theorem 2 of ^1], that the geodesic flow on the ellipsoid with 
distinct semi-axes is strongly nondegenerate, is not true. 

All critical points corresponding to Fq = and F3 = are non-degenerate since Ki{y, y) 
cannot vanish. The type of the points (i.e. whether the eigenvalues are elliptic, hyperbolic, 
or focus-focus) is elliptic on Fq = and F3 = 0, while it changes from elliptic to hyperbolic 
at the degenerate points on Fi = and F2 = 0. This is where the straight line Fi = in 
figure H changes from being the outer boundary of the image of the momentum map to a 
Une interior to the image; similarly for the line F2 = 0. 

The critical points with double roots ri = r2 = d are non-degenerate when they occur 
outside the umbilic curves. In that case the Jacobian matrix of the flow generated by the 
one degree of freedom system + Q('^2(02)) is non-degenerate at A2(02) = d since the 
2nd derivative of Q simply gives 2d, which is non-zero, and gives elliptic type. 

The six corank two points given by the intersection of the lines Fj = and Fj = 
are non-degenerate. The above statements can now be specialised to the plane Xi = 
Xj = Vi = Vj = 0. This is a one degree of freedom geodesic flow on the ellipse. The 
points on this orbit can be easily parametrized and then the eigenvalues become = 
8/tai/(— (ttfe — ai){ai — ai)) where k,l are the other two indices distinct from These 
eigenvalues are real or pure imaginary. The linear combination aDXp^ + jSDXp. has 
eigenvalues aXi and f3Xj, which are different assuming 7^ aj. The combinations that 
occur are elliptic-elliptic (indices 03, 01, 23), elliptic-hyperbolic (02, 13), and hyperbolic- 
hyperbolic (12). 

The topology of the bifurcation diagram is always that of flgure H as long as all the 
semi-axes of the ellipsoid are distinct. As a flrst attempt to understand the degenerate 
cases in flgure Elthe bifurcation diagram is presented for four cases in which the semi-axes 
nearly coincide. From top left to bottom right the cases are close to equal smallest axes, 
equal middle axes, equal largest axes, and equal smallest and largest axes, also denoted by 
211, 121, 112, and 22. It appears as if in the 121 and 22 case the image of the momentum 
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map has only a single chamber, and in the 121-case the image is a triangle. We will 
see that this is not quite the correct answer. From the bifurcation diagrams it is clear, 
however, that the image of the symmetric subspace(s) which correspond to a number of 
collapsing lines needs to be analysed anew, while the results outside this preimage can be 
taken over. 



3 Ellipsoid with equal middle axes 

Consider the geodesic flow on a three ellipsoid with equal middle axes ai = a2- The 
Casimirs Ci and C2 and the resulting Dirac bracket and the Hamiltonian are the same 
as before. However, the integrals Fi and F2 are not defined any more, but the singular 
terms cancel in the sum G = Fi + F2. The other integrals Fq and F3 remain the same. 
The system is invariant under rotations in the {xi,X2) plane and its cotangent lift, i.e. 
simultaneous rotation in the (1/1,2/2) plane. This S0{2) group action is 

^x,y;9) = iS:,y) (10) 

where 

X = (xo, Xi cos 6 — X2 sin 6, xi sin 6 + X2 cos 6, x^) ^^^^ 
y = iVo, yi cos 6-y2 sin 9, yi sin 9 + y2 cos 9, y^) . 

The group action <l> is the flow generated by the angular momentum J = Xiy2 — X2yi, 
which is a global action variable since it generates the periodic flow 

Theorem 3.1. Liouville Integrability The Geodesic flow on the ellipsoid with equal middle 
axes is Liouville integrable. Constants of motion are the energy H = ^ (?/q + yf + y2 + yi), 
the angular momentum J = Xiy2 — X2yi, and the third integral G = Fi + F2 

^_ 2 , 2 , (^0^1 - ^lyof , (^0^2 - X2yo)^ , (3:1^3 - X3yif , {x2y3 - x^yif , . 

Or — +1/2+ + + + l^i/J 

Proof. As in the generic case 2H = Fq + G + F^. So G and H commute because the Fi 
commute in the generic case. In the limit a2 ai the constant of motion (ai — 0:2)-^! 
becomes J^. Hence commutes with H and G, and therefore also J . 

The constants of motion are not independent on the symplectic leaves of the Dirac 
bracket. Their relation is 

Fn G F-i . , 

— + j + — = 0, 13 

ttO CH ttf Ct3 

which is a straightforward limit of the generic relation ^ Fj/aj = using Fi/ ai + F2/ a2 = 
G/a2 - Fi{ai - a2)/(aia2)- 

H, J, G and the Casimirs Ci, G2 are functionally independent almost everywhere on the 
level set Gi = G2 = 0: They are polynomial and independent e.g. at a; = (y^, 0,0,0), 
2/= (0,1,0,0). □ 
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The group action $ has the invariants 



■Ki= + X2, = 2/1 + 2/2, TTs = Xiyi + X2y2: T^i = Xiy2 - X2yi , (14) 

related by 7ri7r2 — 7r| — 7r| = 0. The remaining variables xq, X3, yo, 7/3 are trivial invariants 
of $. The fixed points of $ have Xi = X2 = yi = y2 = 0. When J = = j 7^ the fixed 
points are not in J~^{j) and the reduction by the 5*0(2) symmetry leads to a smooth 
reduced system on J~^{j)/S0{2): 

Lemma 3.2. A set of reduced coordinates (■Coi Ci) ^O; ^i; ^2) is defined on the reduced 
phase space Pj = J~^{j)/S0{2) by the formulae 

Co = 2:0, = V^, 6 = 2:3, vo = yo, vi = V2 = ys- 



The reduced coordinates satisfy the Dirac bracket in R^[C,r7], i.e. 

iik, 6} = 0, {^jk, 77;} = 6ki - — , r]i\ = . 

Dauai JJcxkai 

The mapping R : M.^[x, y] ^ M^[^, rj] is a Poisson map and the reduced system has reduced 
Hamiltonian 

and additional integral 

G^nl+ (^^^0 ~ ^0^^)' + (6r?2 - 6^1)' eg ^ 

«! — ao ai — as Ci V 0(i — ckq cki — 

Proof. Define a set of coordinates on M^[^,r7] as shown. The Poisson property of the 
map R, i.e. {f o R, g o R}^ = {f,g}Q o R follows from direct computation of the 
basic brackets, e.g. {6, 6}6 = + a^j, a^sjs = 0, {6,'72}6 = + ^hy^js = 

—{xlx3 + xlx^) / {Dala3^/xl+~x^) — —$,i$,2/{Daia3), etc. The reduced bracket {, }e has 
the Casimirs Ci = ^o/ao + ^f/ai +^i/as — 1 and C2 = Co^o/«o + + ^2^2/012- The 

relation between the invariants becomes ^\tt2 — rjiCf ~ = and elimination of 7r2 from 
the Hamiltonian leads to the above result. Similarly the integral G can first be written in 
terms of invariants TTj, i = 1, 2, 3, and then the elimination of TTj in addition using tti = 
and — rji^i gives the result. □ 

The reduced system is the "geodesic fiow" on the 2-dimensional ellipsoid with semi-axes 
^/ao, ^/al, ^/a3 and an additional effective potential j^/2^f. By definition ^1 > 0, so that 
the reduced system for \j\ > has only the open half of the ellipsoid as configuration 
space. Since \ j\ > the plane ^ = is dynamically not accessible because — xi — 
X2 — ^ j — 0. Clearly the coordinates Xo,Xs,yo,y3 can serve as local coordinates on 
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the half-ellipsoid, see the singular reduction below. Alternatively ellipsoidal coordinates 
on the ^-ellipsoid can be used to separate the variables. 

To this end define a (singular) coordinate system on M*[xj,?/i] by 

xi = ^icos6', X2 = ^isin6' (15) 

where 9 is the angle of rotation corresponding to the 5*0(2) symmetry group action $. 
The 2-ellipsoid embedded in is defined by Ci = 0. Coordinates (Ao,Ai,A2) are then 
chosen as a confocal ellipsoidal coordinate system [T2] in M^[^], these being the roots z of 



C2 (^2 (^2 

+ + = 1. (16) 

ao — z a\ — z as — z 

Constant Aq defines an ellipsoid, constant Ai a one-sheeted hyperboloid and constant A2 
a two sheeted hyperboloid, where Aq < ao < Ai < ai < A2 < 03. Fixing Aq = gives a 
set of generalised coordinates (Ai, A2, 6) on the three ellipsoid with the middle two semi- 
axes equal. The conjugate momenta are denoted by {pi,P2,Pe), where pe is the angular 
momentum J found earlier. 

Lemma 3.3. The Hamiltonian for the geodesic flow on the ellipsoid with equal middle 
axes in local symplectic coordinates (Ai, A2, ^,Pi,P2;P6») reads 

^ _ 2(ao - Ai)(ai - Ai)(a3 - Ai) ^^ _ 2(ao - A2)(ai - A2)(a3 - A2) ^2 
Ai(A2 — Ai) ^ A2(Ai — A2) ^ 

("0 - ai)(a3 - "1) 2 
2ai(Ai - ai)(A2 - ai) ^' 

The constants of motion are pg and Gi 

~ _ 2(^0 - \i){ai - - Aj) 2 , , («o - «i)(«3 - Q^i) 2 

~ A^ ' ~ 2a,{\-a,) 

where i = 1,2. The integrals G and Gi are related by 

~ ~ _ (ai - a3)(ai - ao)^ a\- a^a-^ ^ 
ij\ + lt2 — la.\ii H 2 Pe- 



a\ a 



Proof. The Hamiltonian in local coordinates is found after performing a cotangent lift 
of the new coordinates, and then expressing the original Hamiltonian in terms of those 
coordinates. The Hamiltonian is separated simply by multiplication with A2 — Ai and 
rearranging to determine Gi and G2- As a result of the separation the momenta pi 
conjugate to A, can be expressed as 



2 ^ _ (17) 

4A(A,) 
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with Q being the analogue of Q in (jHl) given as 



2h{ai-zY+ (18) 

(a3 - a;i)(a;i - ao) . , ^ ■2\ , aoOs-af.^, . 

-{g{ai -z) - J ) + 2 ^ ("1 - ^) 



Oil a 



The relation between the constant of motion G and the separation constants Gi is derived 
by substituting the expressions for the original coordinates in terms of the new local 
coordinates into G, rearranging and expressing in partial fractions. □ 

An analogue of the relation between the constants of motion for the generic 3-dimensional 
ellipsoid (|7j) is given by 

Fo ^ Fs ^ G ^ J' Q{z) 



z — ao z — z — ai {z — ai)^ A{z) 

The separating coordinate system is singular whenever A, equals ak, hence whenever 
= 0. When smooth elliptic coordinates (pi are introduced the singularity ai xi = X2 = 
remains, while those at Xq = and 0:3 = disappear. The umbilical points on the reduced 
ellipsoid {xl,xl) = (ao(«i ~ Q^o)) 0^3(0^3 ~ ch))/(q;3 ~ cto) (coresponding to Ai = A2 = ai) 
are contained in the plane ^1 = 0. 

Lemma 13.31 can be read as singular coordinates: for the full as for the reduced system. 



4 Singular Reduction 

The singular reduction for j = leads to a reduced system on a non-smooth manifold. 
To understand its singularity let us consider as an aside the simple example of the 5*0 (2) 
action $ on the cotangent bundle T*M^ with coordinates {xi,X2,yi,y2)- We can define 
the angular momentum, invariants and reduction map in exactly the same way as we did 
for the case of geodesic flow on the ellipsoid. The reduction map n gives a reduced phase 
space given by the surface 7ri7r2 — 7r| = embedded in M'^[7ri, 7r2, vra]. This surface is a 
cone when j = and smooth otherwise. Considering the inequalities vti > 0, 7r2 > the 
reduced phase space for j = is half of a cone. An alternative description of this reduced 
phase space is obtained by first restricting to any invariant subspace of the Xi-a;2-plane, 
e.g. X2 = y2 = 0. The 5*0 (2) action $ has a residual Z2 action on this plane since 
$(7r)(a;i, yi) = (— Xi, —yi). Therefore the singular reduced phase space {711112 = 7r|} C 
can also be viewed as M.'^[xi,yi]/1>2. This is e.g. the half-plane xi > with the boundary 
yi = identified with itself by {xi,0) ~ (— Xi,0), which again gives a cone. The fixed 
point of the residual Z2 action $(7r) is the origin xi = yi = and it is the singular point 
of the reduced phase space. Yet another representation of the same reduced phase space 
is given by classical polar coordinates xi = rcos9,X2 = rsin6', so that the reduced space 
is the half plane r,pr with r > and pr = {xiyi + X2y2)/r. From the above we see that 
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for j = the reduced space, however, is not a half-plane, but a cone, because of the 
identification pr ~ —pr along the line r = 0. 

The cone 7ri7r2 = 7r| can be diagonalised as a quadratic form hj ni = u + v , 7i2 = u — v , 
and TTs = w. Then the cone is parametrised hj u = r, v = rsin0, and w = rcos0. In 
complex notation z = rexpi(j) the reduction map tti = x^, 112 = y"^, tts = xy can then be 
written as z = i{x — iy)"^ /2. Thus the mapping from M.'^[x,y] \ (0,0) to the cone without 
tip is a double cover. Therefore again the cone is equal to R^/Z2. Moreover the Poisson- 
structure on the cone given by {vti, 712)3 = 47r3, {tti, 713)3 = 27ri, {712,713)3 = — 27r2, is 
mapped into the symplectic structure {711,712)3 = {x^, 7/^)2 = 4:xy = Ans, similarly for 
the other brackets. Moreover, the reduction map is invariant under the Z2 symmetry 
action, and thus the symplectic structure on the plane passes down to the cone. A similar 
argument is valid in the case of geodesic flow on the ellipsoid: 

Lemma 4.1. The singular reduced phase space of the geodesic flow on the 3- ellipsoid 
with equal middle axes and vanishing angular momentum j = is the phase space 
of the geodesic flow on the 2-ellipsoid reduced by the Z2 action S{^o,^i,^2,VojVijV2) = 
(^Oj ~^i5 ^2, ''70, —Vij V^)- Thus it is the geodesic flow on the 2-ellipsoid with a hard billiard 
wall inserted in the ^1 = 0-plane. 

Proof. The S0{2) group action $ does not act freely and so we have to use singular 
reduction to analyse the reduced phase spaces. To do this we use invariant theory. The 
Casimirs for the system, expressed in terms of the invariants (I14|l . are 

Xn VTi , XoVo 7r3 x^y^ „ 

^ + — + ^ = 1, + — + = 0. (19) 

tto Oil 0^3 ao Oil as 

Note that these equations are linear in the invariants. The Hamiltonian may be expressed 
in terms of the invariants as 

H = l {yl + ^2 + yl) . (20) 

The reduced phase space is a subset of M''[xo, yo, X3, y^, ni, 1x2-, T^-i\- It is defined by the two 
Casimirs ()19|) . the relation between the invariants 7ri7r2 — 'n\ = and the inequalities 
7ri > 0, 7r2 > 0. It carries the induced Poisson-structure. To describe this subset we first 
of all eliminate 7ri and 7r3 using to get 



aiU-^l-^l\.2-a\{^-^^^-^\ =f. (21) 
V ao as/ V "0 "3 / 

This is a single equation in M^[xo, 0:3, ?/o; Z/s, 7r2] which defines a four dimensional object. 
Equating the gradient to zero, we find that the equation defines a smooth four dimensional 
manifold except when j = 0. This is the reduced phase space Pj = J~^{j) / S0{2). When 
j = there is a singularity for 7ri = 7r2 = 7r3 = 0. The singular points are given by the 
phase space of the geodesic flow on the ellipse in the 03-plane, which is a cylinder. Thus 
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the singular set of the reduced phase space is itself a symplectic manifold. This symplectic 
manifold is invariant under the flow of the reduced equations, but it is not fixed under it. 

The geodesic flow on the 2-ellipsoid reduced by the Z2 action S is the billiard. Take as a 
fundamental region the half-ellipsoid with ^1 > 0. Then the action of S on the boundary 
^1 = simply flips the sign of rji, which can be viewed as the reflection on the plane 
^1 = with the rule "angle of incidence equals angle of reflection". Moreover, points 
^1 = i]i = ^ fixed under 5* and correspond to orbits that are sliding along (or in) the 
billiard boundary ,^1 = 0. 

To establish the correspondence between the reduced space and the billiard consider 
the slice X2 = 1/2 = through full phase space. This is a geodesic subflow, which is the 
geodesic flow on the 2-ellipsoid with semi-axes ^/aQ, i/aT, v^«3, as described in lemma 
Any motion with j = can be reduced to a motion in this plane by some rotation $(6') 
with constant 9. So locally the flow on this 2-ellipsoid is the reduced system. Globally, 
however, we still have the residual Z2 action S, given by $(7r), to reduce by. 

The reduced bracket in is the original Dirac bracket between xo,x 3,1/0, ys with the 
additional non-zero brackets with 112 given by {i/i, ^2} = 2xiTT2/ As in the trivial 
example above we now show that the mapping from the 2-ellipsoid minus a cylinder to 
the reduced phase space minus the singular set is a Poisson map. When the singular set is 
removed, ^21\ can be solved for 712- After elimination of tt2 only the original Dirac bracket 
between XQ,X3,yo,y3 remains. For the description of the geodesic flow on the 2-ellipsoid 
we use the the variables ^0,^15^2 and momenta rjo,rji,rj2 which satisfy the Dirac bracket 
(0), without implying that they are obtained by regular reduction as in lemma The 
mapping from the symplectic submanifold of M^[^0) ^i? ^2, ^70) Vii V2\ obtained by fixing the 
Casimirs to M^'fxo, X3, y^, 2/3] (without any Casimirs) is simply the projection xq = X3 = 
^2,yo = VOiVs = V2, which preserves the Dirac bracket. However, this mapping is 2 : 1 
since from the Casimirs only can be recovered, but not its sign. 

The meaning of this construction is very simple. Because of the inequality tti > the 
variables xq, X3 are restricted to the interior of the ellipse tti = 0. These variables are good 
local coordinates on the reduced phase space after the singularity (at tti = 0) is removed. 
From the point of view of the 2-ellipsoid this amounts to choosing local coordinates in 
configuration space as the projection of the point onto the ^1 = plane. □ 

As a result the regular reduction described in lemma 13.21 gives the right description 
even in the singular case, when the additional discrete quotient by the Z2 action S is 
included in the picture. In the regular case j ^ the two halves of the ellipsoid that are 
identified by S are dynamically disconnected, because of the effective potential j^/2^^. 
However, when j = the Z2 action is less trivial because its fixed set is now accessible to 
the dynamics, and this fact is crucial in order to obtain the correct singular fibres in the 
next section. 
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5 The Liouville foliation 



We now wish to investigate the topology of the invariant level sets obtained by fixing the 
constants of motion. The energy momentum map is SA4 = {H, J,G) : M ^M.^. Since H 
for a geodesic fiow is homogeneous in the momenta we can fix the energy to, say, h. 

Theorem 5.1. The image of the energy momentum map EM. for constant energy H = h 
is the region in bounded by the quadratic curves (see figure\^ 

9 = h -J , g = h -J . (22) 

ai — as ai{ai — as) ai — ao ai{ai — ao) 

Singular values of the energy momentum map are the boundary curves (elliptic), their 
intersections (elliptic- elliptic), and an isolated singularity at the origin {j,g) = (0,0) of 
focus- focus type. 

Proof. As in the generic case critical points can occur on sub-ellipsoids. On xq = yo = 
the integral Fq = and V-Fq = 0, similarly for = = and F^. In both cases the 
corresponding sub-ellipsoids are ellipsoids of revolution. The image of the critical points 
with xo = = is found using the relation (fT!^ to ehminate F3 in 2H = Fq + G + F3, 
which gives 

2H = G-aJ — -^ 

and hence the first curve of critical values. A similar computation for critical points with 
= ys = gives the other curve. 

These points are non-degenerate because the Jacobian of the fiow generated by Fq 
restricted to the critical points = yo = is given by 0- Evaluating this on the point 
Xq = Xi = X3 = 0, yo = y2 = shows that the eigenvalues of this matrix never vanish and 
are of elliptic type, and hence these critical points are non-degenerate. Similar arguments 
apply to the points x^ = ys = 0. The two corank two points given by the intersection of 
the two curves are also non-degenerate, because the non-zero 2x2 blocks of the Jacobians 
are distinct, so that fiDXp^ + vDXp.^ spans the Cartan subalgebra; the 4 eigenvalues (for 
any point on the critical circles given by Xq = x^, = Vq = y^ = 0) are ±2i^^J2aQh / (ai — ao) 
and ±2ivyj2ash/ (as — ai). This orbit is a relative equilibrium, i.e. a circle in the X1-X2 
plane. The eigenvalues of DXp^^ and DXp^ are elliptic, so at their intersection an orbit of 
elliptic-elliptic type is found. 

Since for Xi = yi = the integrals Fi = and also its gradient vanishes, G = Fi -\- F2 
and its gradient clearly vanishes when Xi = X2 = yi = y2 = 0. Considering the Casimirs 
the solutions set of xi = X2 = yi = 1/2 = is a geodesic fiow on the ellipse in the xq-x^ 
plane. Fixing the energy two critical circles are obtained. On these critical points also 
J = so that the origin in the image (J, G) = (0, 0) is a critical value. 

Moreover VJ = as well, and the Jacobian of Xj has eigenvalues ±i, since its fiow $ 
is a rotation. Finally ^DXg + vDXj spans the Cartan subalgebra; the 4 eigenvalues (for 
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any point on the critical circles) are ±fi^/8aih/{ai — ao){as — a;i)±iz/. These eigenvalues 
form a complex quadruplet, and hence the isolated critical point at the origin is of focus- 
focus type. 

This establishes the existence, non-degeneracy, and type of all critical points. The 
bifurcation diagram is shown in figure El with an isolated singularity at the origin {j,g) = 
(0, 0). The remaining part of the proof shows that there are no other critical points. 

First of all the ellipsoidal coordinates from lemma are used to establish that almost 
all other points are non-singular. These coordinates are non-singular outside any sub- 
ellipsoid = 0. To find critical points in the region of phase space with C, coordinates 
such that all 7^ it is enough to compute the rank of the matrix D{Gi,G2,Pe)^ see 
lemma. 13.31 Since the variables are separated this implies = and h2ai{\s — aiY = 
— («! — ao)(<^3 " OLijPe which is impossible. Thus critical points are contained in the 
coordinate singularities = 0. 

It remains to check the pre-images in full phase space of the sub-ellipsoids = where 
the ellipsoidal coordinates are not defined. When r^j = in addition to = then 
the point is critical, see above. Thus we need to show that all points with = but 
r^j 7^ are non-singular. The ellipsoid xq = is a totally geodesic submanifold, i.e. when 
?/o = every orbit stays inside xq = 0. Conversely, when ?/o 7^ the orbit must leave 
the sub-ellipsoid Xq = 0. Similar for X3 = 0. In general an orbit with x in some sub- 
ellipsoid(s) but y not tangent to these sub-ellipsoids will leave them, and therefore will 
have all Xfc 7^ 0. But there ellipsoidal coordinates are regular, and therefore the original 
point is non-critical, since geodesic motion preserves non-criticality. When ^1 = and 
hence xi = X2 = the condition r/i = is always satisfied by definition, but it does not 
specify yi and ?/2- But the previous argument applies again: if (1/1,1/2) 7^ (0,0) then the 
geodesic will leave Xi = X2 = 0, and therefore the original point is not critical. 

□ 

The bifurcation diagram figure El can be considered as the square root of the diagram in 
figure 121 top right. The line J = is where the lines Fi = and F2 = collapse. However, 
in the limit coming from the generic case the whole line would appear to be critical, since 
it is on the boundary of the image of the energy-momentum map. Recall that the limit of 
^2(0:2 — Oil) (or Fi(q;i — 0:2)) equals J^, but not J. Obviously is singular when J = 0, 
but J itself is not. Thus the square root of figure El top right gives figure El Upon this 
transition most of the critical points along the lines Fi = F2 = disappear, except for 
the isolated critical point, which can be thought of as the remainder of the intersection 
of the two lines and the critical curve corresponding to the umbillic line Ai = A2 (in the 
generic system). Moreover the two corank 2 points at the corners of the triangle in figure El 
top right only have corank 1 after passing from to the "better" constant of motion J. 
Finally, the multiplicity of the regular changes from 2 to 1 for every regular point in 
the image. 

The fibre of a regular value in the image of EM. is a by the Liouville- Arnold theorem. 
We now wish to find the fibres of the energy momentum map at the singular values in 
the bifurcation diagram, in particular at the isolated critical value. 
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Figure 3: Bifurcation Diagram for the ellipsoid with axis ao < ai = a2 < as- 

Theorem 5.2. The singular fibres over the boundary curves of the image of the energy 
momentum map at constant energy, with the exception of their intersections, are two tori 
T^. At each intersection point of the boundary curves the singular fibre is S^. The singular 
fibre over the isolated singularity at the origin is the direct product of and a doubly 
pinched torus T^. 

Proof. At the boundary of the image all singularities are of elliptic type, and hence the 
singular fibre is T'^"'' where r is the corank of the singularity; r = 1 on the upper and 
lower curve and r = 2 at their intersection points. 

The upper boundary with Fq = consists of all orbits in the geodesic flow on the 
ellipsoid of revolution defined by xq = yo = 0. Reduction maps each of this system 
to a relative periodic orbit. The isloated periodic orbit in the 12-plane of the geodesic 
flow on the ellipsoid of revolution corresponds to the extremal points with J = ±v^2ai^. 
Reduction maps this realtive equilibrium to the fixed point ^ = (0, ^/al, 0) on the middle- 
axis of the reduced ellipsoid. A similar statement holds for the lower boundary = 0. 
The isolated periodic orbit in the 12 plane is common to both ellipsoids of revolution. 

The singular fibre for the isolated critical value at the origin of the bifurcation diagram 
contains two circles of critical points in the 03-plane, see above. The critical points are 
non-degenerate and of focus-focus type. From the general theory jHj it follows that it is 
an almost direct product of a pinched torus multiplied by 5*^. Since there are two circles 
of critical points in the singular fibre the number of pinches is 2. 

The fibre over the isolated singularity is complicated as it is not of elliptic type and thus 
contains the non-critical points of the separatrix in addition to the two circles of singular 
points. By lemma EH] the reduced system for j = is the geodesic flow on the 2-ellipsoid 
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Figure 4: The intersection of the preimage of the isolated singular point (j, g) = (0, 0) 
with the Poincare section xs = in reduced phase space. The separatrix is of type C2 
before quotient by Z2. 

quotient by the Z2 action S. Ignoring the quotient the reduced singular fibre consists of 
the unstable isolated periodic orbits in the plane ^1 = and their separatrix. In [2j it was 
shown that the topology of this singular fibre is 6*2x6'^, where C2 stands for two circles 
intersecting in two points. This can be seen from the Poincare section .^2 = 0. Since 
^2 = ''?2 = is an invariant subflow the boundary of the section with 772 > is an invariant 
set and it is the only place where the flow is not transverse to the section. In configuration 
space the section condition is the ellipse in the 01-plane, and it can be parametrised by an 
angle by (^o^i) = (^A^ocos^, sin0). The momentum conjugate to then gives 
the momenta as {riQ,rii) = {y/oosmcp,— y/oi cos (f))ptf,/d where c? = sin^(0) + ai cos^(0). 
The reduced Hamiltonian can be solve for rj2 on the section and thus the integral G can 
be written as a function of (0, p^) on the section: 

2h = -g— + -| + - + cos20 . (23) 

The singular fibre is (7 = (j = was already used in the singular reduction) , and it defines 
two circles winding around the {(f>,p^) cyclinder intersecting at the points (0, 0) and (vr, 0), 
which are critical points of g. This is the "atom" C2 Now the quotient with respect 
to S has to be performed. In the new coordinates the action of S is (0,P(/,) — > (—0, — P</>), 
which fixes ^0 and tjq, but flips the sign of ^1 and rji. This action has two fixed points (0, 0), 
and (tt, 0), so that the fixed points of S coincide with the critical points of g. Reduction 
of the cylinder by the Z2 action gives the "canoe" [S], with two singular points. The 
two singular points are connected by two half-circles. This is C2/Z2, where Z2 acts by 
reflection such that the intersection points are fixed. 

Since the reduced flow is transverse to the section on the singular fibre the complete 
reduced singular fibre is (C2/Z2) x S^. The singular fibre in full phase space is found by 
letting $ act on the preimage of this set under the reduction map. Since the singular 
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circles are fixed under $ they will remain singular circles, while every other point will 
be multiplied by S^. Exchanging the order of the operations, first acing with $ on the 
preimage of C2/Z2 gives a double pinched torus, which is then multiplied by 5*^. This 
action also has a generator, which is a second global smooth action, see below. 

□ 

Following the approach suggested by Nguyen Tien Zung in ^3], one can reformulate the 
last statement of Theorem 15. 21 bv saying that the singularity corresponding to the isolated 
singular point is the direct product of the standard 4-dimensional focus-focus singularity 
with 2 pinches and a non-singular system with 1 degree of freedom. This is a kind 
of "almost direct product" decomposition which can be found for any non-degenerate 
singularity (see [IS]). The fact that in our case the product is "direct" seems to be a 
general property of focus type singularities (Nguyen Tien Zung, private communication). 



6 Actions and Monodromy 

We found that the equations for geodesic flow are Liouville- Arnold integrable, and so the 
fibre over a regular point is a T^. Let Ci, C2, C3 be a basis of cycles on this torus. Due 
to the fact that the variables can be separated as in lemma EISl natural cycles are 

Ci : dXi = d\2 = 0,02-. de = d\2 = 0,C3:de = dXi = 0. (24) 

The adjective 'natural' is used in the technical sense of ^H!, and simply means to consider 
the obvious choice given by coordinate lines of the separating coordinate systems. Natural 
though this may be, it will turn out that the corresponding actions are only continuous, 
but not smooth. The corresponding natural action variables are given by 

= 77- / PedG = Pe, l2 = ^ (f PidXi, I3 = ^ (f P2d\2. (25) 



27r 27r Jc^ l-n jc.. 

The first action is just the angular momentum. The momentum pi in the second and 
third actions is given in lemma ESI together with the polynomial Q in ()18|1 . The integer 
multipliers m2 = 2 and = 2 arise due to the way in which the ellipsoidal coordinates 
are set up over the eight octants in M^, see p. When I2 (respectively I3) is evaluated on 
the upper (respectively lower) boundary of the bifurcation diagram (Figure EI) the natural 
actions for the geodesic flow on the ellipsoid of revolution are found, see appendix. 

The polynomial Q can be factored as Q{z) = z{z — ri){z — r2) where, in order to have 
real momenta pi, p2, 

c^o ri < ai < r2 < as . (26) 
The integrals ()25p are calculated on the hyperelliptic curve given by 



= -^m^. (27) 

[z - aiY 
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Figure 5: Real part of curve w'^ = —z{z — aQ){z — ri){z — r2){z — a^) showing cycles. 

The genus of this curve is one less than for a non-degenerate ellipsoid because the pole 
in A can be divided out. This is plotted in figure El and the cycles can be seen. Note 
that the part of the curve in the negative z range does not correspond to any real motion. 
Writing out the actions in full we have proved the following theorem 

Lemma 6.1. The actions of the geodesic flow on the three dimensional ellipsoid with 
equal middle axes are given by Ii = pe, 

^ ^ ^ dz and h = -l — — ^ . dz . (28) 



TT Jc^ 2{ai - z)w TT 2{z - ai)w 

Note that the constants of motion pg, h and g are implicit in the definition of Q{z) 
given in ()18|1 . The integrand has a simple pole at ai and branch points at ao, as, and at 
Ti and r2 for 7^ 0. The integrals are hyper-elliptic of genus 2 and third kind. The three 
natural actions Ji, I2 and I3, are functions of pe, g and h. However, we will show that I2 
and /s are not differentiable at pe = 0. The derivative of these actions with respect to pg 
is given by 



■ dz (29) 



dlj _ {as - ai)(ai - <yo)pe 
dpe 2-Kai fc, ^i^ti 

Figure ini indicates the poles, branch points and integration paths. Note that as pg — > 0, 
then we find that one of the branch points tends to the pole at ai. 

g < ^ lim r2 = cti, g > ^ lim ri = ai. (30) 

We follow the approach of Cushman jH] and Dullin ^Hj and deform the integration path 
so that the integral may be split up into three or two separate integrals depending upon 
which case is being considered. The integral around the branch points Ci is expanded into 
loops Bi around the poles and then the contributions from the poles R2 are subtracted. 



(31) 

'C2 JB2 JRi JC-i JB3 JR2 
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Figure 6: Complex plane C(s) and choice of branch cuts (a), integration paths for pe for 
ri — > tti (b) and r2 ai (c). (d) and (e) show decomposition of C2 for case (b) and C3 
for case (c). 



Now evaluating the residue of the integrand at the simple pole, we have 

Res r dz = — r-T^ ^- — r (32) 

z=ai [z — ai)w 2[a3 — ai)[ai — aQji\pe\ 
We then have for (7 < 0, 

lim = 0, lim ^-^ = —1, lim ^-^ = 1 (33) 

Pe^o pq pg^o+ dpe pe^o- ape 

and for > 0, 

lim 7—^ = —1, lim = 1, lim — - = 0. (34) 

pe^o+ dpe pe^o- ape pe-^o pe 

In other words we have 

, dli , ^ dli 

hm -7:— t lim 1 = 2,3 

pe^o- dpe po^o+ dpe 

Since I2 and J3 are even functions of pe, differentiability at would imply that the deriva- 
tive is zero at pe = 0. But there is a discontinuity in the derivative here and so the 
natural actions are continuous but not differentiable at pg = 0. By changing the basis of 
cycles locally smooth actions can be found, but they are then globally multi-valued. We 
describe this in the following manner. 
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Let the natural actions for positive pe be represented by /+ = (/i,/2,/3)*, and those 
for negative pe by Ji is odd and J2, /a are even, hence 

I^{-pe) = SI+{pe) (35) 

where S = diag(— 1, 1, 1). We then define unimodular matrices Mj such that J+ and Mj/_ 
join smoothly at pg = 0. This smoothes the actions locally. For continuity on the line 
Pe = it is necessary that below the isolated critical point 

/+ = Mi/„ = M1SI+ = g>Q (36) 

and above 

/+ = M2I- = M2SI+ = M2/+, ^ < (37) 

Therefore (0,/2,/3)* is an eigenvector of Mj with eigenvalue 1 The eigenvector equa- 
tions show that M,- must have the form 



Mi= \ 1 \ (38) 




and since Mj G SL{3, Z) as it is unimodular, we must have 6i = 1. 

Now note that ^ = and |^ = — sgn(p5i) for (7 > and so we may find Mi because here 

M.^ = 1^ (39) 

ope ope 

which implies that Ki = —2 and Pi = 0. 

For g < we have |^ = and |^ = —sgn{pg) and by 

M.§^ = 1^ (40) 

ope ope 

we find that K2 = and (32 = —2. 

The monodromy matrix M for a counter clockwise cycle around the origin [pg, g) = [0, 0) 
is given by M = {M2S)-^{MiS), which gives 



M = 2 10. (41) 




21 



By a final unimodular change of basis, defined by TMT ^ = N, where an appropriate 
choice of T is 



we put the monodromy matrix into normal form and have proved the following theorem: 

Theorem 6.2. Monodromy The obstruction to the existence of smooth global action vari- 
ables for the geodesic flow on the ellipsoid with equal middle axes is monodromy. The 
monodromy matrix has normal form 



This result is consistent with the general theory of non-degenerate singularities of in- 
tegrable systems developed in [TT], [T3], According to Nguyen Tien Zung [T3], each 
non-degenerate singularity can topologically be presented as an almost direct product 
of "basic" singularities. In our case, this is just the direct product of the 4-dimensional 
focus-focus singularity with two pinched points and a non-singular system with one degree 
of freedom (see the last statement of theorem 15. 2j) . It is easily seen that the monodromy 

A 2\ 

matrix in such a situation is decomposed into two blocks: ( ^ ^ j which correspond to 

the focus-focus singularity with two pinches (see pi], ^Hj) and the trivial 1-dimensional 
block. Up to a change of basis, this is exactly the matrix from theorem 16. 21 

A The geodesic flow on 2-enipsoids 

Here we briefly describe the well known classical situation of the 2-ellipsoid embedded 
in M^. There are three critical periodic orbits obtained by intersecting the ellipsoid with 
any coordinate plane Xi = 0, i = 0,1,2. The middle plane Xi = contains the four 
umbilic points defined by Ai = A2 = ai where Aj are (algebraic) elliptic coordinates on 
the ellipsoid. Introducing elliptic functions Aj(0j) gives a double covering by a torus with 
coordinates (pi branched over the umbilic points. To analyse any motion that hits the 
umbilic points these ellipsoidal covering coordinates cannot be used. It turns out that the 
only motion that ever crosses the umbilic points are the two unstable periodic orbits in the 
xi = plane and their separatrices. The image of the momentum map for fixed energy is a 
line segment with three corank 1 non-degenerate critical values. The endpoints correspond 
to the stable orbits in the xq = and X2 = plane. There is another critical value in 
the middle of the interval corresponding to the orbits in the plane Xi = 0. The fibre of 
this point contains the two unstable orbits connected by a heteroclinic separatrix. The 
topology is that of two circles intersecting in two points multiplied by a circle. This can 




(42) 
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be seen by considering the Poincare section xq = with xq = i/q > on the section. The 
boundary of the section Xq = Uq = arc itself geodesies, otherwise the flow is transverse 
to the section. The topology of the section with > is that of a flnite cylinder, i.e. 
an annulus. The restriction of the second integral G to the surface of section and to 
constant energy gives a function with two minima (corresponding to the two geodesies in 
the X2-plane) and two saddles (corresponding to the two geodesies in the Xi-plane). The 
separatrices intersect the section along two curves wrapping around the cylinder once and 
intersecting in two points, like = ± cos 0. 

If two semi-axes are the same we obtain a prolate or oblate ellipsoid of revolution. For 
this system it is an elementary exercise to compute the non-trivial action. For the oblate 
case where the longer axes are equal, ao < ai — a2 the action is given by 



1 

2^ 



Psds (43) 

C2 



where 



Pl=(2h- "° ^^ A i<^l + (o.^-aoy) _ 

The essential integral is a function of the ratio p — ao/ai and the scaled angular momen- 
tum j — j/ \j2ha\ only. Expressing li in Legendre normal form gives 

\-^^ = um-'jj^{^i\^\ u'-i-ni-p), (45) 

where S and IT are Legendre's complete elliptic integrals of the second and third kind 
respectively and the modulus and parameter are 

e^i-pu-^, p^^e/i-p. (46) 

For the oblate case p > 1 and thus k'^ < 0. Upon replacing U by iU the above formula 
for the action also holds in this case. 



References 

[1] A. V. Bolsinov and A. T. Fomenko. Orbital classification of geodesic fiows of two- 
dimensional ellipsoids. The Jacobi problem is orbitally equivalent to the integrable 
Euler case in the dynamics of a rigid body. Funktsional. Anal, i Prilozhen., 29(3) :1- 
15, 96, 1995. 

[2] A. V. Bolsinov and A. T. Fomenko. Integrable Hamiltonian systems. Chapman & 
Hall/CRC, Boca Raton, FL, 2004. Geometry, topology, classification. Translated 
from the 1999 Russian original. 



23 



ll II 

2.5r 2.5r 




Figure 7: Actions for ellipsoids of revolution: (i) ao < «i = oti-, (ii) a\ = a2 < as, where 
ao = 1, ai = 2, as = 4, /i = 1. 
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